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Abstract 

Calculations of cross-sections for the gluon-photon production of and B* mesons 
as well as for the production in ep-interactions at HERA are performed in the leading 
order of perturbative QCD. We show that the nonabelian structure of QCD leads to an 
essential forward-backward asymmetry in these processes. 

1 Introduction 

A significant progress in the understanding of the mechanisms for the production of heavy 
quarkonium consisting of two heavy quarks with different flavors, for instance, (6c)-system, 
has been achieved during past three years 0, ^, |^ . One clarified that for the production 
in e"'"e~-annihilation, the differential cross-section has the simple factorized form in the limit 
of high energies {M%^/ s ^1)0, 

dzdpT dkx 

where z = lE^^j \fs^ o^i is the Born cross-section for the 66-pair production, is the trans- 
verse momentum of 6-quark, is the transverse momentum of B^ meson, and Di^b^{z) is 
interpreted as the fragmentation function of 6 —> i?c + A. In the language of Feynmann 
diagrams, the latter fact means that in the high energy limit one can conveniently choose 
a gauge, namely the axial gauge, with the unit n vector directed along the four-momentum 
of the 6-quark, so that the contribution of the diagrams where the cc-pair is produced off 
the fe-quark line is dominant. Using this fact, one analytically calculated the fragmentation 
functions for the bound states in S-, P- and D-waves 0. 

The production of (6c)-quarkonium in hadronic, photonic and photon-gluon interactions 
occurs in a more complex way and it can not be described by eq. (1) in the entire kinematical 
region even at high energies M^^/s <S 1. The number of diagrams for the mentioned processes 
is greater than that in e"'"e~-annihilation (one has 36 diagrams for the gluonic interactions, 
20 for the photonic collisions and 24 for the photon-gluon subprocess). So, one can find no 
gauge to restrict the number of dominant diagrams by those of the fragmentation kinds , i.e., 
by diagrams like 15 or 18 in Fig. 1, where the tree process g'-f —>■ B^. + X is presented. The 
origin of the problem is that the recombination diagrams also contribute to the cross- sect ion, 
for example, diagram 5 in Fig. 1 is of this kind. This contribution is essential up to the 
transverse momentum values equal to 35 GeV. At pt > 35 GeV, the regime of fragmentation 
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(1) becomes dominant, as one must expect by the theorem on the hard process factorization 
0. However, the latter fact underhnes the calculations of the contributions, determined by 
the complete set of diagrams. The complete perturbative calculations can lead to a possible 
interference of various contributions in some kinematical regions. So, detailed analysis for the 
gluonic production of Be mesons at high energies was performed by us ||^. We found that the 
effect of destructive interference between the fragmentation diagrams and the ones containing 
the three-gluon vertex is important, so that in the region of z close to 1 and at moderate 
Pt values the complete set of diagrams leads to smaller cross-section in comparison with the 
one calculated according to (1) |^. We also found that for the photon-photon interaction 
the interference is absent. These investigations lead us to the conclusion that destructive 
interference takes place if the interaction possesses the nonabelian character. Therefore, it is 
interesting to study the production of Be meson in the gluon-photon interactions, wherein the 
mentioned interference would result in the forward-backward asymmetry, which is described 
in this paper. 

Unfortunately, a straightforward comparison of experimental measurements with the frag- 
mentation picture is problematic since the condition M^^/ s ^ 1 is generally not satisfied in 
partonic subprocesses, and the dominant contribution to the cross-section is integrated over 
the region close to the threshold of the meson production in the subprocess. In this work 
the cross-section of Be production in the HERA experiments will be presented. Despite the 
dominance of the kinematical region near the subprocess threshold, the effect caused by the 
presence of the three-gluon vertices, i.e. the nonabelian texture of QCD, turns out to be 
significant. 



2 Calculation technique 

The amplitude for the B^ meson production, A^'^^'', can be expressed through the amplitude 
of four free quarks production, T^'^^{pi, k{q)), and the orbital wave function of the B^ meson, 
\E''^'^(q), in the meson rest frame as 

A''^' = J T'''{P^, Mq)) ■ (^"'^(q))* ■ (2) 

where J and jz are the total spin of the meson and its projection on z axis in the B^. rest frame, 
respectively; L and Iz are the orbital momentum and its projection; S and Sz are the sum of 
quark spins and its projection; C^^f^ are the Clebsch-Gordan coefficients; pi are four-momenta 
of Be, b and c, q is the three-momentum of b quark in the Be meson rest frame; k{q) is 
the four-momentum, obtained from the four-momentum (0, q) by the Lorentz transformation 
from the Be rest frame to the system, where the calculation of T'^*^(pj, A;(q)) is performed. 
Then, the four-momenta of b and c quarks, forming the Be meson, will be determined by the 
following formulae with the accuracy up to |qp terms 

Pe = ^PB.-Hci), 

where rrib and me are the quark masses, M = rrib + rrie, and Pb^ is the Be momentum. Let us 
note that for the S-wave states it is enough to take into account only terms independent of q. 
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The 24 diagrams, contributing to the gluon-photon production of Be are shown in Fig. 1. 
The product of spinors f jMc, corresponding to the b and c quarks in the amphtude T^^"" (pi, k{q)) 
of eq.(S), should be substituted by the projection operator 



v{T) = Vm ^ r — ^ , (4) 



where F = 7^ for S* = 0, or F = i*{PB^,Sz) for S = 1, where e{PB^,Sz) is the polarization 
vector for the spin-triplet state. 

The fact that the 6c-system is produced in the colour- singlet state simplifies the texture of 
colour matrix. The latter has the only nonzero eigenvalue for the given process, so that one 
can significantly reduce the time for the computation of the matrix element. 

The ^-wave production amplitude can be written down as 

where Rs{0) is the radial wave function at the origin, so that ^^(O) = R{0)/\^ATr and 

Rs{0) = 

and the quantity /s^ is related to the leptonic decay constants of pseudoscalar and vector Be 
states 

{0\UO)\V) = ifvMve,, 
(0|J5^(0)|P) = ifpp,, 

where and Jz^{x) are the vector and axial- vector currents of the constituent quarks. In 

the lowest order in one has / = /v, /p- The account for hard gluon corrections in the first 
order in a, ISI, ISH results in 




f = h 
f = h 



f m2 — nil ^2 8 
1 m 

TT \m2 + nil ITT'l 3 

^ _ af /m2 - mi ni2 



(6) 
(7) 



71 \ni2 + nil ^1 

where mi^2 are the masses of the quarks composing the quarkonium. For the vector currents 
of quarks with equal masses, the BLM procedure of scale fixing in the "running" coupling 
constant of QCD gives (see ref. |TT|) 



The estimates of the Jb^ value within the potential models and QCD sum rules are somewhat 
uncertain, Jb^ = 460 ± 100 MeV p. Recent estimates within the QCD sum rules and in 
the framework of lattice computations give /^^ = 385 ± 25 MeV and fs^ = 395(2) MeV, 
respectively. Taking into account radiative corrections, the value of Jb^ is about 470 MeV 
|T^, which will be used in the following estimates. Since the cross-section is proportional to 
/^^, the corresponding value of the cross-section for a different value of J'b^ can be obtained 
by rescaling with the multiplicative factor {fBc/fs 
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3 Results of calculations 



In the leading order of perturbation theory, the total cross-section for the photon-gluon pro- 
duction of Be is proportional to a ~ """3 |^(0)P, where |\I/(0)|, and masses of quarks are 

b c 

free parameters of the model. The calculations are performed with the following values of the 
parameters 



a 


= 1/133, 


as 


= 0.24, 


TUb 


= 4.8 GeV, 


TUc 


= 1.5 GeV, 




= 470 MeV. 



The results for the gluon-photon production of S-wave levels in the (6c)-quarkonium [Be and 
B*) are presented in Fig. 2 and Tab. 1. For the sake of comparison, the predictions of 
the fragmentation mechanism are also shown in Fig. 2. One can see that the fragmentation 
mechanism underestimates the cross-section, much the same way as in the photonic and gluonic 
interactions producing Be (or B*) 0, ^. 

The distributions da/dz for both the Be meson moving in the gluon hemisphere (i.e. at 
Pz > 0, where pz is the projection of the momentum on the beam axis in the c.m.s. of the 
colliding particles) and Be moving in the photon hemisphere, are shown in Fig. 3. According 
to the fragmentation model, the mesons are produced symmetrically in both the hemispheres, 
so that the distributions over z are described by the expression D{z) ■ a^i/^ aX Pz > Q as well 
as at < 0. However, one can see that the real picture of production is more complex. First, 
the large enhancement compared to the predictions of the fragmentation model is found in 
the region of low z for both the hemispheres. Second, a significant asymmetry takes place 



for z > 0.5. Both effects are caused by the following reasons. As was shown in |T^, the first 
reason is the contribution of diagrams like the diagram 20 in Fig. 1. In the leading order of 
perturbation theory the latter contribution corresponds to the "resolved" photon split into the 
pair of charmed quarks, so that one of those recombines with the produced 6-quark to form 
the _Bc-nieson. Taking into the ratioof the quark charges {Qe/QbY = 4, this mechanism causes 
the increase of -Bc-production cross-section in the photon hemisphere. The second reason is 
the destructive interference between the 6-quark fragmentation and the contribution caused by 
the three-gluon vertices appearing in the splitting of initial gluon inside the cone of the latter 
(see diagrams 13, 14, 17, 18, 22 and 3, 4). This is why the essential decrease of i?c-production 
cross-section occurs in the gluon hemisphere. 

The distributions da/dpx for the forward and backward hemispheres {pz > and Pz < 0, 
respectively) are presented in Fig. 4 for the Be state. One can see that these distributions 
differ from the ones obtained in the fragmentation model at low transverse momenta. In 
the region of high pt, where the fragmentation mechanism dominates, the deviation becomes 
small, as one expects. 

Comparing Figs. 3 and 4, one can draw the conclusion that the difference between the 
production cross- sect ions in the two hemispheres is most significant at large z and low pt- Our 
calculations show that the same conclusion is valid for the production of the excited vector 
state B*. 

The considered interactions between the monochromatic beams of gluons and photons do 
not occur in a real experiment. To clarify whether one can observe the asymmetry in ep or jp 
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experiments or not, the process of ep-interaction at HERA {^/s = 314 GeV) is considered. For 
this purpose, the cross-section of the gluon-photon subprocess has been convoluted with the 
structure functions of both the photon in electron and the gluon in proton, respectively. The 
equivalent photon parameterization in has been chosen as the photonic structure function. 
The CTEQ4 functions have been used for the gluon distributions ||15|| . 

The rapidity distribution normalized to unity is shown in Fig. 5 for the production of 
B*. For comparison, we have also presented the distribution with a constant matrix element, 
|Ap = const. One can see that the differential cross-section for B* production is strongly 
displaced in comparison with the distribution calculated with the constant matrix element. 
Therefore, the effect of destructive interference, which has been found at high energies of the 
(77 — s> Bc{B*) + X subprocess, survives in the real conditions of interaction at HERA, too. 
On the other hand, in the fragmentation regime at high pt the asymmetry must tend to zero. 
The latter practically occurs in Fig. 6, wherein the cross-section da/dy at pt > 30 GeV is 
shown in comparison with the distribution for the constant matrix element at the same cut 
off, Pt > 30 GeV. 

Summarizing the results above, one can conclude that the effect of the forward-backward 
asymmetry in the i?^*)-meson distribution da/dy is basically related with the nonabelian 
texture of QCD, and this asymmetry is determined by the three-gluon vertex. 

As has been mentioned above, the total cross-section of Be and B* production at the 
HERA facilities is about to 2.5 pb including charge-conjugate states. The vector B* state 
decays to the ground Be meson through the emission of a 7-quantum in the Ml-radiative 
transition, which has the branching ratio close to 100% 0. So, at the HERA luminosity 
C = 16 ■ 10^°s~^cm~^ the calculated value of the cross-section will give only 400 events with 
the Be meson per year. The increase of total energy by three times, which correspond to the 
combination of LEP and LHC accelerators, will result in the enhancement of the asymmetry 
as well as the total cross-section, which is estimated to have a value of 6 nb at ^/s = 660 GeV 
(see Fig. 7). At the luminosity 10^^s~^cm~^ one can expect about 6 ■ 10'^ i?^*-* events per year. 

Taking into account low branching fractions for the exclusive decay modes of B^ (for 
example, Br{B^ — >• J/ipe^v) ~ 2.7%) and the efficiency for the detecting of decay process 
(e ~ 0.1), one can see that the expected number of events can turn out to be insufficient for 
the investigation of the described asymmetry in the ep-production of B^ meson. 

The study of these nonabelian effects is a more real problem for the production of Dg 
meson. As we can see above, the asymmetry is caused by quite similar deviation from the 
fragmentation mechanism due to the pro cess- dependent contributions, which have a more 
rapid decrease with the growth of the transverse momentum of the produced meson. For the 
i?^*^ mesons these sub-leading terms can be reliably calculated in perturbative QCD. In the Ds 
production, higher orders in as can be significant. Nevertheless, we believe that the sub- leading 
non-fragmentational contributions have the same feature as in the B^ production: the yield in 
the photon hemisphere is enhanced and the one in the gluon hemisphere is suppressed due to 
the nonabelian interactions of gluons. One has to note that the quark-mass ratio determining 
the asymmetry value in the model under consideration is one and the same for the Ds and 
Be mesons, and hence the effect of destructive interference can be experimentally observed at 
HERA, where the energy scale is given by the ratio of M^Jy/s, which is approximately equal 
to Mb^I \fs at the energy of LEP and LHC combination. The distribution for the Ds-v[veson 
rapidity at the HERA energy is shown in Fig. 8. 

A more detailed study of the interference effect requires the investigation of the contribu- 
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tions by individual diagrams in a proper gauge. The latter consideration will be presented 
elsewhere. 
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Fig. 1. The diagrams contributing to the gluon-photon production of (6c)-quarkonium 
in the leading order of perturbative theory. The c-, 6-quarks, gluons and photons are 
denoted by the thin, thick, spiral and wavy lines, respectively. 




Fig. 2. The dependence of the gluon-photon cross-section for the production of B* (sohd 
triangle) and Be (empty triangle) in comparison with the predictions of the fragmentation 
mechanism (solid and dashed curves, respectively). 
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Fig. 3. The da^^ jdz for the meson moving in the gluon hemisphere (dashed histogram) 
and for the one moving in the photon hemisphere (dotted histogram) at the energy of 
gluon-photon interaction equal to 100 GeV. For the sake of comparison, the prediction of 
the fragmentation model is also presented (solid curve). 
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Fig. 4. The differential cross-section for the gluon-photon production of Be versus the 
transverse momentum at the energy of 100 GeV. The notations are the same as in Fig. 
3. 
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Fig. 5. The differential cross-section normalized to unit area over the rapidity of pro- 
duced in ep-interactions at ^/s = 314 GeV (solid histogram). For the sake of comparison, 
the normalized distribution with the constant matrix element is shown by the dashed 
histogram. The rapidity y is defined in the c.m.s. of the colliding beams. 
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Fig. 7. The same as in Fig. 5, but with the energy of ep-inter actions, E — 660 GeV. 
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Fig. 8. The differential cross-section normalized to unit over the rapidity of D* produced 

in ej9-interactions at ^/s = 320 GeV (solid histogram). For the sake of comparison, 
the normalized distribution with the constant matrix element is shown by the dashed 
histogram. 



